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We investigate the effects of the new scalars in a two-Higgs-doublet model on 
the weak magnetic dipole moments of the fermions at the Z peak. Proportionality 
of the Yukawa couplings to the fermion masses, and to tan (3, makes such effects 
more important for the third family, and potentially relevant. For the r lepton, 
the new diagrams are suppressed by v T = 2sin 2 0^—1/2, or by powers oim T /Mz, 
but may still be comparable to the SM electroweak contributions. In contrast, 
we find that the new contributions for the bottom quark may be much larger 
than the SM electroweak contributions. These new effects may even compete 
with the gluonic contribution, if the extra scalars are light and tan (3 is large. 
We also comment on the problem of the gauge dependence of the vertex, arising 
when the Z is off mass shell. We compute the contributions from the new scalars 
to the magnetic dipole moments for top-quark production at the NLC, and for 
bottom and r production at LEP2. In the case of the top, we find that the SM 
electroweak and gluonic contributions to the Ztt vertex are comparable. The 
new contributions may be of the same order of magnitude as the standard-model 
ones, but not much larger. 
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1 Introduction 



Thus far, the nature of the symmetry-breaking sector of the standard model (SM) 
remains largely untested, with a fundamental scalar yet to be found. In fact, the 
number of Higgs multiplets is not predicted and must be determined by experiment, 
just like the number of fermion families. Moreover, the success of supersymmetry 
coupling-constant unification suggests that one might have several Higgs doublets, 
since the minimal supersymmetric standard model (MSSM) requires two doublets. In 
this article we consider a two-Higgs-doublet model (THDM), which has physical scalars 
H ± , A , H° and h°. 

While these particles remain undiscovered, we may look into the virtual effects that 
they induce in several phenomena. For instance, from the neutral-meson mixings one 
can set limits on the flavour- changing scalar vertices. These limits are so stringent that 
Glashow and Weinberg JTJ and, independently, Paschos ||, introduced the concept of 
Natural Flavour Conservation, implemented through a discrete symmetry. We will 
assume that one doublet couples to the right-handed up-type quarks, while the other 
doublet couples to the right-handed down-type and charged-lepton fields. This is the 
so-called model II, which includes the scalar sector of the MSSM as a particular case. 
In model II, the Yukawa couplings are proportional to the fermion masses, so that one 
should look to the third family for noticeable effects. 

Another interesting source of information comes from the radiative corrections to 
the gauge couplings of the fermions. On the one hand, the Z couples to a fermion / 
through a vector (vf) and an axial- vector (a/) coupling: 



«(P-) \f ( v f((l 2 ) ~ a M 2 )lh)] v(p. 



2s\yCw 

These couplings occur already at tree level in the SM, with 

v f = T 3f -2Q f s 2 w , 

a f = T 3/ . (2) 

Here, sw and cw are the Weinberg angle's sine and cosine, and Qf and T 3 f are the 
fermion's charge and third component of the weak isospin, respectively. The mo- 
mentum of the Z is q = p_ + p + . Therefore, precise measurements are required to 
disentangle the loop effects. On the other hand, the anomalous weak magnetic dipole 
moments (WMDM) fif, defined to be the couplings of the Z to / of the form 

IP r i 

u( P -)\»f(q 2 )icT^qv\v( P+ ), (3) 



2m/ 

arise only at loop level, making them preferred tools in the search for physics beyond 
the SM. 

In this article, we compute the WMDMs of the r lepton and of the bottom quark in 
the model described above. We separate the SM contributions from the ones involving 
new scalar particles, and study the conditions under which the latter are numerically 
important. In particular, we reproduce the SM computation of the r-lepton's WMDM 
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in Ref. ||. We give both analytic formulas for the WMDMs, and also their numerical 
values. The analytic formulas can also be applied to compute the usual magnetic dipole 
moments of the fermion with the photonF]. 

The article is organized as follows. In section 2, we discuss the WMDM of the r 
lepton in the THDM. In section 3 we present the results for the bottom quark. Section 
4 is devoted to a discussion of the problem of the gauge non-invariance of the WMDM 
when the Z is off mass shell. This is crucial for the top quark, but also for the other 
fermions, when they are produced in colliders such as LEP2 or the NLC We give some 
numerical values relevant for these colliders. We draw our conclusions in section 5. An 
appendix contains the analytic expressions for the Z and 7 magnetic dipole moments 
(MDMs), induced at one-loop level in model II. 



2 The WMDM of the r 

One expects the r lepton to be a leading candidate in the search for new physics 
through WMDMs. This is due to the absence of gluon contributions, and to the fact 
that the energy and angular distribution of the decay products of the r can be used to 
extract information on the spin density matrix of the r pairs. Several groups have used 
this method to isolate the dispersive and absorptive parts of the weak electric dipole 
moment (WEDM) g % and of the WMDM g |, [| . An analysis, including all the form 
factors in terms of the spin density matrix of the r, is given in Ref. H. In particular, 
the transverse polarization of the r (within the collision plane) measures the real part 
of the WMDM and the imaginary part of the WEDM, while the normal polarization 
(perpendicular to the collision plane) measures the imaginary part of the WMDM and 
the real part of the WEDM. These are gauge- invariant, observable quantities, as long 
as the external particles are on mass shell: p 2 _ = p 2 + = m 2 , and q 2 = M§. 

In the SM, as well as in model II, the WEDM is multiloop suppressed (since it must 
involve the quark mixing matrix), and therefore these polarizations measure directly 
the WMDM. The calculation of the WMDM of the r in the SM has been performed 
in Ref. ||, who found 

/z^ M (M§) = (—2.10 — i 0.61) x 10~ 6 , (5) 

which is well below the expected experimental limit of 10~ 4 . 

A typical diagram contributing to the magnetic dipole moments is shown in Fig. 1. 
We will denote the different diagrams by the particles running in the loop, starting 
with the particle A in between the external fermions, and proceeding counterclockwise 
to particles B and C. We denote by the SM Higgs scalar, and by a (x) the charged 
(neutral) Goldstone bosons. 

1 In our notation, the tree- level coupling of the photon to / is ieQfj^, while the magnetic dipole 
moment is defined as in Eq. (||) . The usual definition for the photon coupling 0] , 



ieQfu{p-) 



7^i(9 2 ) + ^^V^(g 2 ) 

2m f 



v(p+) , (4) 



yields the following relation with our photon magnetic dipole moment: (if — Q/i*2- 
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The various contributions are suppressed by different powers of m T . One power is 
always there due to the definition of fi T [See Eq. (|3])]. Extra powers come from the 
Yukawa couplings (in the case of the scalars and of the Goldstone bosons) and, in the 
diagrams in which the Z couples in the loop to two particles with the same spin, from 
the mass insertion required to reproduce the chirality structure of the WMDM. The 
latter suppression is avoided when the B and C particles in Fig. 1 are a scalar and a 
vector (or vice- versa), that is, in the T<pZ and uaW diagrams. However, in the SM as 
in model II, the suppression thus avoided is offset by the proportionality of the Yukawa 
couplings to the fermion masses. That is not the case in the most general THDM, 
without Natural Flavour Conservation. 

We have performed a complete computation of fi T in the SM, carefully checking 
its gauge invariance. Each of the four sets of diagrams in Fig. 2 is separately gauge 
invariant when the external particles are on mass shell. (One must also use the tree- 
level relation cw = Mw/Mz to explicitly verify gauge invariance.) This is also the case 
for the diagrams involving the SM Higgs particle (the ones in Figs. 3a and 3b, with 
h — > <f>). The expressions for the MDM in the 't Hooft-Feynman gauge are listed in 
the appendix. We agree with the results in Ref. |§, except for a few signs, and for the 
small T(pZ and rZip contributions, which do not affect significantly the final numerical 
results. 

For the numerical computations, we used the following values. At q 2 = Mf , a em = 
1/128 and s^y = 0.232. (Notice that, due to the large cancellation in v T = 2$^ — 1/2, 
some of the numerical results are very sensitive to the input value of sfy.) We use 
M z = 91 GeV, M w = c w M z , and m T = 1.777 GeV. We take the neutrino to be 
massless. Using this input, we get the values for the contributions to the r WMDM 
displayed in Table 1. We have grouped the diagrams into the gauge- invariant classes 



ABC 


..ABC 

H'T 


ry TT 


(3.1853 — i 1.2713) x 1(T 7 


Ztt + x TT 


(4.1316 + zl. 9133) x 1(T 8 


Wvv + ovv 


(-9.8734-i5.0748) x 10" 7 


vWW + vaW + vWa + voo 


-1.4733 x 10" 6 


Total 


-2.1008 x 10"^ - % 6.1548 x 10" 7 



Table 1: Standard- model contributions to the WMDM of the r 



shown in Fig. 2. Our final result agrees with that of Ref. 0. 

We stress that all the results presented in this article are based on exact formulas, 
which we have evaluated using two completely different programs, as a cross-check. In 
many cases, the numerical results may be guessed at from the analytical expressions 
in the appendix by taking the relevant integrals to be dominated by the largest mass 
scale in the diagram. In the appendix, we show explicitly that this 'guesstimate' of the 
integrals is correct for the jtt diagram. Experience shows that it also works for many 
other diagrams, making it a valuable tool in understanding the relative magnitudes of 
the different contributions. 

We now consider the scalar-particle contributions, both in the SM and in the 
THDM. The relevant diagrams are the ones shown in Fig. 3, and their contributions 
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are given in the appendix. There are four types of scalar contributions: those with a 
neutral scalar and a Z (Fig. 3a), those with a neutral scalar (Fig. 3b), those with a 
neutral pseudo-scalar (Fig. 3c), and those with a charged scalar (Figs. 3d and 3e). In 
the SM, only the first two types of diagrams are present. They are functions of the 
Higgs mass which we display in Figs. 4 and 5, for the ripZ + rZ0, and 0rr diagrams, 
respectively. We have taken all the scalar masses to be greater than 50 GeV. One can 
see that the SM Higgs contributions are never larger than about 2% of the non-Higgs 
contributions. 

In model II, the contribution of diagram 3a to the WMDM is proportional to 
cos 2 a + tan (3 sin a cos a, for H°, and to sin 2 a — tan (3 sin a cos a, for h°, where a is 
the mixing angle in the CP-even mass matrix, and tan (3 is the ratio between the two 
vev's. In the following we shall take 1 < tan/3 < 70 and let a take any value. The 
proportionality factor is the function plotted in Fig. 4. We find that this contribution 
may increase or decrease the real part of the WMDM by about 80%, for a scalar of 
mass 50 GeV. Of course, this diagram does not originate an imaginary part since, at 
q 2 = M§, no cut can be made. 

Similarly, the amplitude generated by diagram 3b is proportional to the function of 
the mass of the scalar shown in Fig. 5. The proportionality factors are (1+tan 2 f3) times 
cos 2 a for H°, or times sin 2 a for h°. The appearance of tan 2 (3 makes this contribution 
potentially large. It turns out that the contribution to the real part of the WMDM is 
negligible, but that the contribution to its imaginary part can make it to about one 
third of the SM value, for a scalar of mass 50 GeV. This happens despite the fact that 
this diagram has an extra m 2 /M| suppression factor. 

The contribution from the pseudo-scalar in diagram 3c is equal to tan 2 (3 times a 
function of the mass of A displayed in Fig. 6. We see that, if that mass is 50 GeV and 
tan 2 /3 is large, this may decrease the imaginary part of the WMDM by about 25%. 
The influence on the real part of the WMDM is negligible. 

Finally we turn to the contribution from the charged scalar involved in diagrams 
3d and 3e. We show in the appendix that such contributions are, in general, linear 
combinations of 1, tan 2 /3, and cot 2 /?. However, due to the zero mass of the neutrino, 
in this case only the term proportional to tan 2 (3 is non- vanishing. The factor of pro- 
portionality is displayed in Fig. 7 as a function of the mass of the charged scalar. We 
see that both the real and the imaginary parts of the charged-scalar contribution may 
make it to about 1 x 10~ 6 , if the mass of that particle is 50 GeV, and tan/3 is maxi- 
mal. This is a contribution of the same order of magnitude as the SM one, but of the 
opposite sign. 

In conclusion, we find that all the new diagrams arising in model II may have an 
impact on the WMDM of the r lepton, if tan (3 ~ 70 and the mass of the scalar particles 
is ~ 50 GeV. They may be comparable to the SM contributions, but never much larger 
than them. This is partly due to the masslessness of the neutrino, which eliminates 
some terms from the general expressions, and partly due to the proportionality to the 
vector coupling of the tau, v T , which is approximately —0.036. This will not be the 
case for the quarks. 
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3 The WMDM of the bottom quark 



The measurement of the WMDM of the bottom quark is more problematic since its 
polarization is affected by hadronization. Still, one may look for its influence on the 
Z — > bb observables. Currently, this yields bounds of the order of 1CT 2 ||. 

Another possibility is brought about by the fact that the initial spin of the bot- 
tom quark is retained in the polarization of the A& baryons that are produced directly 
fT0| , P| , although there might be substantial depolarization induced by the A b 's pro- 
duced from the decay of heavier b-baryons f[0|, [12| . This polarization can be studied 
by measuring the energy spectra of the charged-leptons or the neutrinos [T3| , 
produced in the semileptonic decays of the A b . 

In the calculation of the WMDM of the bottom quark, we have used the values 
of M z , M w , sw and a em given in the previous section. In addition, we take m b = 5 
GeV, m t = 174 GeV, and as(M z ) = 0.117. The values obtained are listed in Table 2. 
Contrary to what happened in the case of the r, for the case of the bottom, all the 



ABC 


..ABC 

n 


gbb 


(3.5764 -i 1.9382) x 10" 4 


^ybb 


(1.9900 -i 1.0785) x 10~ 6 


Zbb + xbb 


(3.0041 + % 1.3393) x 10~ 6 


Wtt + att 


-2.2064 x lO" 6 


tWW + taW + tWa + too 


-6.4791 x 10~ 6 


Total 


(3.5394 -% 1.9356) x 10" 4 



Table 2: Standard-model contributions to the WMDM of the bottom 



gauge-invariant electroweak contributions are of the same order of magnitude, except 
for the Higgs contributions, which amount to less than 1% of the overall electroweak 
result. However, the WMDM of the bottom quark is dominated by the gluon diagram, 
shown in Fig. 8, which is about two orders of magnitude larger than the electroweak 
ones. 

Still, the new scalars in model II may be important since the two factors limiting 
their impact on the WMDM of the r lepton are not present in the case of the bottom 
quark. In fact, the Yukawa-coupling factors in the vertices are now often factors of 
enhancement instead of suppression, because of the large coupling of the top, and 
also of the bottom in the large-tan (3 limit. In addition, the vector coupling v b = 
2/3 sin 2 6 W - 1/2 « -0.345 is an order of magnitude larger than v T . 

The shape of the contributions from diagrams 3a and 3b as functions of the scalar 
masses are similar to those of Figs. 4 and 5, although the numerical values are different. 
We find that the diagrams in Fig. 3a may increase or decrease the real part of the 
WMDM by some 35%, if one of the scalars of the THDM is very light while the other 
one is very heavy, and if tan j3 is very large. Similarly, the h°bb or H°bb diagrams in 
Fig. 3b may be important for the imaginary part of the WMDM, increasing its SM 
value by about 70% if tan/5 is large and the scalars are light. 

As before, the pseudo-scalar contribution is proportional to tan 2 (3 and, if tan 2 {3 = 
5000 and m^o = 50 GeV, that contribution is maximal, reaching (—3 + 9i) x 10~ 5 . 
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Finally, we look at the contributions from the charged scalar H ± . These are real 
since we are assuming that the mass of the charged scalar is larger than Mz/2. As 
we have mentioned above, these contributions are linear combinations of 1, tan 2 /3, 
and cot 2 /?. We find that these contributions are, at best, two orders of magnitude 
smaller than the QCD contribution, for tan 2 f3 = 5000 or tan 2 f3 = 1 (in which case it 
is dominated by the term proportional to 1). 

So, contrary to what happened for the r lepton, for the bottom quark only the 
neutral-scalar and pseudoscalar contributions may change the SM value for the WMDM 
appreciably. One should note the remarkable fact that these new scalar contributions 
may compete with the strong gluonic corrections. 



4 Magnetic dipole moments at q 2 > M 2 



z 



The CDF and DO Collaborations |15| have shown that the top quark is quite heavy. 
Because of this fact, problems arise in the definition of an appropriate WMDM for 
the top, since in order to produce a top pair the Z must be off mass shell. The root 
of the problem lies in the gauge dependence of any form factor arising in the Lorentz 
decomposition of the Zff and 7// vertices, when the gauge boson {Z or 7) is off 
mass shell ||16|| . Indeed, the MDMs are unequivocably defined (gauge- invariant and 
observable) only when the incoming gauge boson is on mass shell (q 2 = for the 7, 
and q 2 = Mf for the Z), i.e., when the physical process is dominated by the single- 
gauge-boson exchange. 

Of course, this problem does not exist solely with the top quark. In particular, 
it occurs whenever one studies the radiative corrections to fermion-pair production at 
colliders with energy beyond the Z mass, such as LEP2 or the NLC. 

Recently, the pinch technique |I7| has been used to construct gauge-independent 



MDMs |]18| . This technique consists in reorganizing the usual Feynman- diagram ex- 
pressions into portions that are manifestly gauge independent. One extracts from the 
box diagrams those gauge- dependent pieces which are kinematically equivalent to the 
7// or to the Zff vertex corrections. Those pieces offset the gauge-noninvariance of 
the vertex corrections. 

One may now calculate these quantities in any gauge as long as the pinch con- 
tributions are consistently identified in the box and vertex graphs. However, in this 
particular case, the computation in the 't Hooft-Feynman gauge is most convenient 
since, in that gauge, the MDMs do not receive any contributions from the box dia- 
grams. Indeed, such contributions from the box diagram arise from the longitudinal 
terms in the gauge-boson propagators, which are absent in this gauge. This means 
that, in the 't Hooft-Feynman gauge, the MDMs only receive contributions from the 
vertex corrections themselves. Thus, the fact that kinematically the top quark produc- 
tion requires q 2 > 4:m 2 , would be circumvented by a judicious definition of the form 
factors when the vector boson is off mass shell. 



However, other authors [[nj have claimed that the improved off-shell vertex func- 
tions are not uniquely defined. The basic problem is that one can still shift gauge- 
independent pieces between the various diagrams. This ambiguity means that the 



7 



gauge-invariant quantities obtained from the pinch technique are not observable by 
themselves. Still, they may be useful in determining which new physics effects may be 
important, and where to look for them. 

Conscious of these problems, we have computed the contributions to the MDMs of 
the top quark with the 7 and with the Z at yfq 1 = 500 GeV (the center-of-mass energy 



expected for the next e + e~ linear collider, NLC) following the prescription of Ref. |l 
We have also computed the MDMs of the bottom and of the the r at yfq 1 = 200 GeV 
(relevant for LEP2). We do this in order to have an idea of the sensitivity of those 
MDMs to the new physics in the THDM. Our calculations are to be considered as 
preliminaries to a more general work in which the box diagrams should also be taken 
into account, in the computation of suitably defined physical observables. 



4.1 The top quark at the NLC 

For the NLC, we have taken as = 0.096, a em = 1/127 and sin 2 9w = 0.240. These 
values are obtained by a SM renormalizat ion-group running of these parameters from 
their measured values at q 2 = M|. In Table 3, we present the SM contributions to 
the WMDM in the Zit vertex at yfq 1 = 500 GeV. Notice that, contrary to the naive 



ABC 


..ABC 
rt 


gtt 


(-2.6170 + ^4.5493) x 10~ 3 


jtt 


-7.1549 x 10~ 5 + i 1.2438 x 10" 4 


ztt + xtt 


-2.4585 x 10~ 4 - % 1.4502 x 10" 3 


Wbb + abb 


-4.2241 x 10~ 4 + % 1.2600 x 10" 3 


bWW + baW + bWa + baa 


-2.2436 x 10~ 3 + % 8.7342 x 10" 4 


Total 


(-5.6004 + i 5.3569) x 10~ 3 



Table 3: Standard-model contributions to the WMDM of the top, 



at y/q 1 = 500 GeV 



expectation, the electroweak sector gives contributions comparable to the gluonic one. 
Qualitatively, this difference with respect to the case of the bottom quark is due to the 
fact that here the largest mass scale that may dominate the integrals is never much 
larger than m 2 , while the mass of the bottom satisfies m 2 . M^,m 2 . 

Now we can discuss the effect induced by the new scalars of the THDM. If the 
new scalar particles can be produced directly at y/q 2 ' = 500 GeV, their effect on the 
MDMs will be to generate new non-vanishing imaginary parts due to unitarity. We 
are more interested in studying a possible scenario, in which they only have virtual 
effects. Therefore, we are looking for the possibility of sizeable effects induced by a 
charged scalar with mass larger than 250 GeV (meaning that no H + H~ event has been 
detected) and neutral scalars bounded by the following constraints: m^o H o > 409 GeV 
(no scalar-Z event); and m h o H o + m^o > 500 GeV (no scalar-pseudoscalar event). 
The formulas for the various contributions are listed in the appendix. Notice that for 
the top quark there is a tan j3 cot (3 interchange with respect to the formulae for a 
fermion with T^f = —1/2. 
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We find that the contribution from diagram 3a, which is real since we have taken the 
scalar masses to be greater than 409 GeV to avoid direct production, is now maximal 
for tan/? ~ 1 but may only reach 10% of the SM value for scalar masses close to 409 
GeV. This value decreases very rapidly as one goes away from 409 GeV into higher 
masses, as a consequence of the threshold effect that exists close to production. 

Similarly, diagrams 3b and 3c take on their maximum values for tan/3 ~ 1. This 
will limit their impact, although the extra m 2 /Mf prefactors enhance these contribu- 
tions. We find that these contributions never go beyond 15% of the SM result, in their 
imaginary and real parts. 

Finally, the charged-scalar contributions permit a test of both regimes: tan (3 ^> 1 
and tan /3 ~ 1. In the first regime the dominant piece is the one proportional to tan /3, 
which may reach two thirds of the SM value in its real part, but only 12% of the SM 
value in its imaginary part. For the tan/5 ~ 1 regime, one may reach around 20% of the 
SM real part and decrease the SM imaginary part by about 15%. These numbers are 
obtained for a charged scalar with mass 250 GeV, but decrease very slowly for higher 
values of the mass. 

The SM contributions to the MDM coupling of the •yit vertex at \fq 1 = 500 GeV 
are listed in Table 4. The gluonic contribution for this vertex is enhanced with respect 



ABC 




gtt 


-8.2790 x 10~ 3 + % 1.4392 x 10~ 2 


jtt 


(-2.2635 + ^3.9349) x 10~ 4 


Ztt + tft 


-6.4295 x 10~ 5 — i 3.5312 x 10" 3 


Wbb + abb 


(-1.4404 + i 4.2711) x 10" 4 


bWW + baW + bWa + baa 


-1.4978 x 10- 3 + % 7.0128 x 10~ 4 


Total 


(-1.0212 + i 1.2383) x 10~ 2 



Table 4: Standard-model contributions to the MDM of the top with the 7, 

at y/q* = 500 GeV 



to the previous one by Qt/xt ~ 3.16. This fact makes it the dominant contribution to 
this MDM. 

Again the contributions from the new neutral scalars may, at most, reach 20% 
of the SM value, being more important for the imaginary part of the MDM. (Recall 
that the photon does not have a contribution from diagram 3a.) In contrast, the 
charged-scalar contributions may reach two thirds of the real part of the SM value, for 
tan 2 /3 = 5000 and a charged-scalar mass near the threshold for H + H production. 
These contributions are not as large for tan 2 /3 ~ 1 where one can only reach 15% of 
the real part of the SM value. 

4.2 The r lepton and the bottom quark at LEP2 

The calculations of the WMDM and the electromagnetic MDM for the bottom and 
t at LEP2 follow those in the previous section. At LEP2, with ^fq 1 = 200 GeV, we 
use as = 0.106, a em = 1/127.5 and sin 2 9 W = 0.236. We remind the reader of the 
ambiguity in defining appropriate MDMs off the intermediate vector boson mass shell, 
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and report our results only qualitatively. We take these as a hint of where to look for 
important virtual effects of the new scalars. We shall take the charged-scalar masses 
to be greater than 100 GeV, together with the constraints m h o H o > 109 GeV and 
mh°,H + wia > 200 GeV. 

Phenomenologically, the r may be more interesting, due to the possibility of a clean 
measurement of its polarization. The SM values are: (—1.1142 — % 3.4060) x 10~ 6 for 
the WMDM, and (1.3347 - % 2.6968) x 10~ 6 for the MDM with the photon. These 
values are comparable with those obtained for the WMDM at Mf . 

Due to the arguments explained in section 2, the contributions from the new scalars 
in the THDM only affect the MDMs in the region of large tan (3. In the case of the 
WMDM, the contribution from diagram 3a may be very large. For a mass of the scalar 
near the threshold of 110 GeV, one gets two times the SM real part. All the other 
contributions are negligible, except for the charged-scalar ones which can reach 25% 
(8%) of the real (imaginary) part of the SM result, but with the opposite sign. This 
occurs for tan 2 /3 = 5000 and charged scalar masses of 100 GeV. 

For the photon vertex all the diagrams may give relevant contributions, for tan 2 f3 = 
5000 and light scalars. Diagram 3b may have an imaginary piece as large as 40% of 
the SM one. The pseudo-scalar correction may reduce the real part of the SM result 
by 60% and its imaginary part by 30%. Finally, the charged scalar may increase the 
real part by 25%. 

For the bottom quark at LEP2, we have found the following results: the standard- 
model WMDM is (7.6143 - i 4.0685) x 10" 5 , while the standard-model MDM with the 
7 is (5.8639 — i 3.4638) x 10~ 5 . These values are about a factor of five smaller than 
the WMDM at q 2 = Mf . Both MDMs are dominated by the QCD correction to the 
vertex. 

In the case of the Z vertex, there is a substantial contribution to the real part coming 
from diagram 3a, which may amount to three times the gluonic one, for a neutral scalar 
close to 109 GeV. On the other hand, the dominant new scalar contributions to the 
imaginary part come from diagram 3b, which may increase the gluonic value by 70%, 
and diagram 3c, which may decrease it by 50%. The pseudo-scalar diagram may 
decrease the real part of the WMDM by 25%. All these results are for tan 2 /? = 5000. 

For the photon vertex, the largest corrections to the imaginary part that one may 
obtain come from diagram 3b (increase of 70%) and diagram 3c (decrease of 50%), for 
large tan (3. Also, the pseudo-scalar diagram may decrease the real part by 25%. 

5 Conclusions 

We have performed a complete analysis of the magnetic dipole moments of the fermions 
in the SM and compared them to the corrections that may be induced by the virtual 
scalars in a THDM with a discrete symmetry, the so-called model II. For the r lepton, 
the masslessness of the neutrino and the smallness of v T suppress the new contributions 
which may, nevertheless, be as large as the SM electroweak corrections. For the bottom 
quark, these suppressions disappear and the new physics yields MDMs much larger 
than the electroweak ones. They may even compete with the gluonic correction. We 
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have also presented the results for the MDMs defined by the pinch technique, which 
are relevant for the top quark produced at the NLC, and for the r lepton and for 
the bottom quark produced at LEP2. This definition is ambiguous, but it is useful 
in identifying candidate situations in which to look for new physics. In so doing we 
find some interesting results. The pure electroweak radiative corrections to the Zti 
vertex at NLC are as large as the gluonic corrections, although that is not the case for 
the jti vertex. The diagrams in Fig. 3a can give contributions to the real part of the 
WMDM of the bottom and r two or three times larger than the SM values. In most 
other cases, precise measurements of the magnetic dipole moments would be required 
to disentangle the new effects. Those effects are maximal for light scalars, allowing for 
tests of the region of parameter space with large tan /3, both in the case of the r and 
of the bottom. 
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A Formulas for the magnetic moments 

In this appendix we give the magnetic dipole moments of a fermion / induced at one- 
loop level in the THDM. We present the results in the 't Hooft-Feynman gauge. They 
hold for any fermion, and are written for both the Z and 7 magnetic moments, for any 
value q 2 of their squared momentum. However, some of the results are gauge-invariant 
only if the gauge bosons are on mass shell, as we have found by computing them in a 
general 't Hooft gauge. 

For simplicity, we introduce the definition 

^ABC = ^J_ h ABC (6) 

Arc 

where A, B and C are the particles in the loop, in the order displayed in Fig. 1. 
Moreover, Xf and yj are defined as 

_ ( v f /(2s w c w ) \ ( a f /{2s w c w ) \ , 

Xf -{ Q f ) ' Vf ~ { ) ' (7) 

where the upper lines hold when the exterior gauge boson is the Z, while the lower 
lines hold for the 7 MDMs. The quantities Vf and Qf are given in Eq. (0). The letter 
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/ refers to the fermion whose MDM is being calculated, and % to its SU(2)^-doublet 
partner. 

In order not to clutter our formulas, we will omit possible mixing-matrix elements, 
notably in the quark sector. The inclusion of those mixing-matrix elements is clearly 
trivial. Furthermore, our formulas are valid for any fermion. In the text, they have 
only been used to compute the MDMs of the fermions of the third family, because those 
are the ones for which the contributions from the extended scalar sector are larger. 

We use for the integrals the notation of Ref. || : 

[loo; I(tu] (m A ,m B ,m c ,q) 

d^k [1; k^] k^k v ] , 



m 2 (k 2 — m 2 A )[[k — pS) 2 — m? B ] [(k + p + ) 2 — ttiq] 
They are decomposed as 

f = ( P _-p + ri 10 +(p_+p + ri n , 

IV = (pW + +p^_)I 21 + (pW_+p^ + )I 22 

+ W + - Vtp v _) Vx + gT ho ■ (9) 



These functions are trivially related to the Passarino-Veltman functions |pO|l . If m B = 
me-, In and I2-1 vanish. Also, notice that, from I^ u , only the combination J 2 i — I22 
appears in the expressions for the MDMs. 

We have obtained the numerical values of these integrals in two independent ways. 
On the one hand, FeynCalc was used to reduce them to functions of the scalar integrals 
A , B and Co, for which there are explicit formulas. The results were checked against 



the ones obtained with the aid of the ff routines of van Oldenborgh [21 . 



A.l SM contributions, without the Higgs scalar 

In this subsection we include the contributions from diagrams involving only SM par- 
ticles, except the Higgs scalar. The latter contributions will be included in the next 
section, where we discuss the scalar sector of the THDM. We find 

= 8x f Q 2 f [I 2 i - I 22 - ho] (0, m/, m f , q) . (10) 

For quarks only, there is a similar contribution with a gluon g instead of the photon: 

tfff = 8x f C F ^[I 2 i- J 22 - I 10 ] (0, m f , m f , q) , (11) 

where Cp = 4/3. For these contributions, the integrals can be computed analytically 
and one obtains the explicit formula 

1 f 1 dy 



[hi - I22 ~ ho] (0, m f , m f , q) = — J 



Aq 2 Jo y(y — 1) + m 



2 /„2 



logg, forg 2 <0; 

(12) 



2q 2 5 \ log + in , for q 2 > 4m 2 f , 
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where 



5 = 



1 - 



Am 
I 2 



f 



(13) 



The region < q 2 < Am 2 is unphysical. In the limit q 2 = one obtains l/(Am 2 ), 
reproducing the Schwinger term. Thus, this integral is dominated by the largest mass 
scale: it is of order 1/m 2 when Am 2 3> \q 2 \, and of order 1/q 2 when q 2 3> Am 2 . This is 
a common feature of these integrals, which is crucial to obtain the correct decoupling 
limits. 

For the diagrams in Fig. 2b we find 



b zff 



foXff — 



b w c w 



Xfm 2 



s 2 M 2 

b W lvl W 



Xfvj (hi - I22 - ho) + x f a 2 f (hi - I22 - 5/io + 2/oo) 
+2y f v f a f (hi - I22 ~ 3/io + ho)] (M z ,m f ,m f , q) , 

[hi - I22] (M z ,m f ,m f ,q) . 



The sum of these two contributions is gauge-invariant. 
Similarly, the sum of the two contributions 



■Wii 



Vi 



[hi ~ I22 ~ 3/io + loo] ( M w, m { , q) , 



1 



, , r , { \m)(xi + yi) + m 2 (xi - yi)] (hi - I22 ~ ho) 

S W M W 



+ 2xim 2 ho] (M w ,mi,m,i,q) , 



is gauge-invariant. 
Finally, 



(14) 



(15) 



iWW 



LicrW 



= s 



iWc 



cot 9 W 
1 

cot (29 w ) 
1 



— [2/21 - 2/ 22 + I 10 ] (m h M w ,M w ,q), 
s w 



1 



s 2 M 2 
s w m w 



m 2 f + mf) (hi - hi ~ ho) 



+m 2 (/oo-2/i 
tan (9 W ) \ 1 



[mi,M w ,M w ,q) , 



2 s 2 



ho(mi,M w ,M w ,q) . 



(16) 



Once again, only the sum of these four contributions is gauge invariant. Moreover, 
gauge invariance in this case only occurs if, either q 2 = M§ and one uses the upper line 
for the coefficients (for the case of a Z as external gauge boson), or q 2 = and one uses 
the lower line (for the case of a 7 as external gauge boson). This is the reason behind 
the gauge noninvariance of the MDMs at arbitrary q 2 . Gauge invariance of a physical 
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observable (such as the e + e~ — > ff cross section) is restored if we add to Eqs. (16) the 
result of the computation of the H^-box diagram . 

In Eqs. (|16|), S is +1 if the fermion / has T 3 f = —1/2, and —1 if it has T 3 f = 1/2. 
This we may write as S = —2 T 3 f. The factor S arises from the anti-symmetry of the 
three-gauge-boson vertices under the W + «-> W~ interchange. 



A. 2 The scalar contributions 



We now present the contributions to the MDMs that involve scalar particles in the 
loop. The letters h, H and A refer to the neutral scalars h°, H° and A of model II, 
respectively. The charged scalars are denoted by the letter C . The corresponding 
Feynman rules may be found, for example, in Ref. f2"2fl . 

The formulas we present are for a fermion with T 3 f = —1/2, such as the bottom 
or the r. For a fermion with T 3 f = 1/2, such as the top quark, one must make the 
following substitutions 

sin (3 <-> cos (3 , sin a <-> cos a. (17) 



We find 



LiCC 




tan/?) 2 + (mi/ tan (3) 2 (I 



'21 



hi ~ I 



10) 



+m 2 (2/ 10 - loo)} K, M c , M c , q) , 



lCU 



b Aff 
b Hff 
b hff 



s 2 M 2 
■>w m w 



Xf-mj 



s w lv± w 

,2 



| (m f tan (3) 2 (xi + yi) + (m^/ tan/3) 2 (xj - ?/j 
-2^m 2 / 10 }(Mc,m i ,m i ,g), 

tan 2 (3 [I 2 i - I22] {M A ,m f ,m f , q) , 



21 



'22 



iio) 



4^4 00^ [/21 - 122 - 2/lo] (M ^' m/ ' g) 



x/m 2 sin 2 a 



sLM 2 cos 2 /? 



[hi - hi ~ 2/io] (M h ,m f , m f , q) . 



These diagrams contribute both to the photon and to the Z vertices. 

In addition, there are a few diagrams contributing exclusively to the Z vertex. In 
some of them the external Z couples either to H° and A , or to h and A , in the 
loop. Those diagrams do not yield any contribution to the WMDM. In other diagrams 
the external Z couples either to a Z and a H°, or to a Z and a h°, in the loop (see 
Fig. 3a). These diagrams have the same chiral properties as the MDM operator and, 
as a consequence, their contributions do not get (rrif/Mw) 2 pre-factors. The results 
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are 

b fZH + b fHZ = Vf^ ^ CQg 2 q; _|_ tan (3 sin a cos a) [I w + In] (m f , M z , M H , q) , 

b fZh + b fhZ = _!V_ (sin 2 a _ tan/5 

sin a cos a) [ii + hi] (m f ,M z ,M h ,q) . (19) 

s w c w 

To obtain the SM Higgs scalar contributions from the expressions in this subsection, 
one deletes all the diagrams with H ± , A , or H°, and sets 

P = a + ir/2 , (20) 

in the diagrams with h°. 
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FIGURE CAPTIONS 



Figure 1: Typical Feynman diagram contributing to the magnetic dipole moments 
at one loop. 

Figure 2: Electroweak one-loop vertex corrections to the MDMs, excluding the 
diagrams with the Higgs scalar. 

Figure 3: One- loop vertex corrections due to the scalars in the THDM. The scalars 
H°, h°, A and H ± are denoted by the letters H, h, A and C, respectively. For the SM, 
only diagrams 3a and 3b exist, having the Higgs particle as the intermediate neutral 
scalar. 

Figure 4: Plot of the contribution of the diagram in Fig. 3a to the WMDM of 
the funtion of the mass of the neutral scalar. This is to be multiplied by 

cos 2 a + tan [3 sin a cos a for H°, and by sin 2 a — tan f3 sin a cos a for h°, in the THDM. 

Figure 5: Plot of the contribution of the diagram in Fig. 3b to the WMDM of the 
r, as a funtion of the mass of the neutral scalar. This is to be multiplied by (1 +tan 2 f3) 
times cos 2 a for H°, and times sin 2 a for h°, in the THDM. 

Figure 6: Plot of the contribution of the diagram in Fig. 3c to the WMDM of the 
r, as a funtion of the mass of the neutral pseudo-scalar. This is to be multiplied by 

tan 2 (3. 

Figure 7: Plot of the sum of the contributions of the diagrams in Figs. 3d and 3e 
to the WMDM of the r, as a funtion of the mass of the charged scalar. This is to be 
multiplied by tan 2 j3. 

Figure 8: Gluonic one-loop vertex correction to the MDMs of the quarks. 
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